A Lattice Model for Colloidal Gels and Glasses 
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We study a lattice model of attractive colloids. It is exactly solvable on sparse random graphs. 
As the pressure and temperature are varied it reproduces many characteristic phenomena of liquids, 
glasses and colloidal systems such as ideal gel formation, liquid-glass phase coexistence, jamming, 
or the reentrance of the glass transition. 
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There is a growing experimental and theoretical inter- 
est in the study of colloidal suspensions of mesoscopic 
particles immersed in a polymeric solution. Their prop- 
erties can be controlled precisely via chemical or physical 
manipulations, and the inter-particle interaction tuned 
from a very short-range depletion-induced attraction to 
a long-range Coulombic repulsion l] . This makes colloids 
important both in the terms of basic scientific research 
and industrial applications 0] . Colloidal suspensions also 
offer new scenarios of dynamically arrested structures, 
which are the main focus of this work. 

At high densities these systems exhibit the same glassy 
dynamics as a system of hard spheres [3]. The lack 
of available space at densities p w 0.58 gives rise to a 
phase called repulsive glass. If a short-range attraction is 
present, ergodicity may also be broken, resulting in the 
so-called attractive glass. When the range of attraction 
is less than « 10% of the sphere diameter, the interplay 
between these two mechanisms can lead to a remarkable 
reentrance of the glass transition line, with the (ener- 
getically driven) attractive glass dominating at low tem- 
perature and the (entropically driven) repulsive glass at 
high temperature. This reentrance was first predicted 
by the mode coupling theory and later confirmed by ex- 
periments and simulations [4J. Interestingly at low tem- 
peratures the attractive glass transition moves to much 
lower densities. Such a low-density arrested state is often 
identified as a colloidal gel 0] - a low-density disordered 
arrested state which does not flow but possess solid-like 
properties such as a yield stress. However, the dynamical 
arrest at low densities and its exact relation to the glass 
transition is still very poorly understood, and our work 
also sheds light on this aspect. 

In the low density - low temperature regime, a sys- 
tem of hard-core particles with a short-range attraction 
undergoes a phase separation between a colloid-rich and 
a colloid-poor phase [6]. The interplay between phase 
separation and gelation is also a widely discussed topic. 
Typically, when quenching from high temperature, the 
phase separation boundary is crossed first. When the 
glass transition is crossed at lower temperatures, the 
spinodal decomposition induces the formation of dense 



arrested regions, a situation called arrested phase separa- 
tion [7J . Another phenomenon can be obtained by reduc- 
ing the phases-coexisting regions to even lower densities 
and temperatures, by decreasing the average coordina- 
tion number of particles, as proposed in [3,0] • Under 
these conditions the phase separation does not intervene 
in a wide region of densities and the system can be equi- 
librated by annealing down to very low temperatures, 
leading to a state of matter called ideal gel [10( • 

Theoretical progress has been hampered by the lack 
of microscopic models able to roughly locate and de- 
scribe the colloidal gel transition. In this letter we dis- 
cuss an exactly solvable yet realistic model characterized 
by valence-limited interactions that generalizes the lat- 
tice glass model introduced by Biroli and Mezard (BM) 
in [ll|]. The model is motivated by the classical short- 
range depletion attraction but also by recent studies of 
a family of patchy Q or limited- valency potentials [l2j ]. 
It allows to reproduce the behavior of a broad variety 
of complex materials that are usually simulated at high 
computer cost. Under different conditions we observed 
in particular the liquid-glass coexistence, the reentrant 
glass transition and the formation of an ideal gel. This 
leads to an unifying description of the structural arrest 
in different amorphous and disordered systems, ranging 
from hard-spheres to colloids and gels. 

The model - We consider a system of particles on a 
lattice. Each site i carries n.$ = or m = 1 particles, with 
a short-range attractive interaction giving a unit energy 
gain to all pairs of nearest-neighbors particles. Motivated 
by the concepts of geometrical frustration and valence- 
limited interactions [H, 0, [TH, we follow BM 11] and 
restrict the occupation by a hard constraint: a particle 
cannot have more than £ occupied neighbors. The parti- 
tion function (where /3 is the inverse temperature and fi 
the chemical potential) reads 

Z(fi,f3)= e"£«^ n4+/3 (1) 

allowed {n} 

where the sum is restricted over all configurations {n} 
satisfying the hard constraints. For (3 = 0, the attraction 



2 



is turned off and the model reduces to the BM model 11 1 . 
In real (continuous) space, the presence of the maximum 
possible number of neighbors I gives rise to an entropic 
loss [lj|. This effect is, however, inhibited on the lat- 
tice and we thus also consider a variant of the model 
where each sites with exactly i occupied neighbors induce 
a penalty e~ Sp in eq. |T]). This entropic repulsion is only 
relevant at large temperatures and is small compared to 
the attractive potential at low temperatures. 

The solution - We shall focus on the exact solution 
of the model on large c— regular random graphs (where 
the number of neighborhooding sites, connectivity, of 
each vertex is exactly c > I). Such graphs are locally 
tree-like, i.e., there are no cycles up a distance scaling like 
logiV from a typical site. The topological disorder and 
the geometrical frustration are introduced by very long 
loops that disfavor the crystalline order and let emerge 
more easily the glassy phases. This corresponds to a 
mean- field like approach to the problem [3] . 

In the low density-high temperature region the Bethe 
approximation gives an exact solution on tree-like lat- 
tices. For s p = (but the generalization is straightfor- 
ward) it is obtained from the fixed point of 
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where ip e /ip s /ip u are the probabilities that a node is 
empty/saturated/unsaturated when one of its neighbor- 
hooding sites is removed from the graph, Zq is a normal- 
ization constant ensuring that ip e + tjj s + tjj u = 1. The 
free energy density J = log Z/N then reads 



J(p,f3) = log 
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The density of occupied sites is given by p — d fl J(p,p). 

As the density increases, or the temperature decreases, 
a dynamical glass transition arises at Pd(P) where the 
free energy landscape divides into exponentially many 
valleys [15[ . Each valley corresponds to a different Gibbs 
state a with a free energy density J a . This can be recog- 
nized from the divergence of the point-to-set correlation 
length, or, cquivalently, by the emergence of a nontrivial 
glass solution described by the one-step replica symmetry 
breaking approach [IE>[lg, II) HI- The complexity S(J) 



(sometime called configurational entropy) is defined as 
S(J) = lagJ\f(J)/N ) where Af(J) is the number of states 
with free energy density J. It can be computed via its 
Legendre transform <I>(to, p, (3) defined by 



e N*(m,n,0) = y-] e % 



NJ a (p,0) = f dJg A'[mJ( AI ,/3)+S(J)] 
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FIG. 1: (color online) High temperature behavior of the model 
with c = 3, I = 1 and no entropic penalty, prcp — p is the 
difference between the packing fraction, or density, and its 
maximum achievable value (random close packing), p is the 
chemical potential. The static equilibrium solution (bold/red) 
undergoes a phase transition at pK = 0.5725 and converges 
to the densest packing prcp = 0.57574 at p — oo. How- 
ever, the dynamics falls out of equilibrium at pd = 0.5708 
and the packing fraction following a slow annealing can be in- 
creased up to at most 0.5731 < p.j < 0.5748. The dotted/blue 
curves are Monte-Carlo annealings for N = 10 5 with rates 
dp = 10~ 4 , 10~ 5 , 10" 6 and 10" 7 (from top to bottom) and 
the (dashed/black) line pDyn is the iso-complexity approxima- 
tion for the infinitely slow annealing (see the text). The inset 
shows the same plot on a logarithmic scale. 

where the sum runs over all the different Gibbs states a 
and the Legendre parameter m is e quiv alent to the Parisi 
parameter in the replica theory [l5l [l9j . The cav- 
ity method derives a set of self-consistent equations for 
3>(m,p,,/3) [H i . T hev can be solved analytically in the 
limit p — > oo ll[. In general, however, these equations 
are solved using the numerical population dynamics tech- 
nique of [15|. The case m = 1 also allows significant sim- 
plification [17|. We also checked self-consistently, using 
the method of [l8[, that the one-step replica symmetry 
breaking solution is stable, and the presented equilibrium 
solution is thus generally considered exact. 

The glass and jamming transitions — We first revisit 
the behavior the of the model in the high temperature 
limit {[3 — 0) where it reduces to the original BM [llT |. 
and consider the case c = 3,£ = lfor concreteness. FIG.Q] 
depicts in bold/red the difference between the largest 
possible density of occupied sites /?rcp and the equilib- 
rium density p as a function of the chemical potential 
p. In the equilibrium solution, the model undergoes a 
Kauzmann transition at px — 0.5725 where the com- 
plexity vanishes and the system becomes a an ideal glass 
(a "static" glass phase) [14J. The density of the glass 
can be in principle further increased in the equilibrium 
solution, and the densest possible packing is reached at 
Prcp = 0.57574 (Random Close Packing, first computed 
in ]n|.) 
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FIG. 2: (color online) Three representative phase diagrams of the model in the temperature (T) — packing fraction (p) plane. 
The black/continuous curves are the colloid-rich/colloid-poor coexistence lines III, and the corresponding spinodal lines are 
black/dashed curves. The orange/dashed-dotted curves correspond to the percolation of physical clusters, and intercept the 
critical point of the phase separation. The dynamical Td(p) (blue/dotted curve) and the Kauzmann Tk{p) (red/continuous 
curve) glass transition lines are also shown. The black filled area is beyond the random close packing fraction prop- Left 
panel: c = 10, I — 3 and no entropic penalty. An ideal gel appears at low enough temperature (T < 0.1) and low volume 
fraction (0.1 < p < 0.4). Central panel: — 6,1 = 4 and no entropic penalty. A phase coexistence between a colloid-poor and 
an attractive colloid-rich glass phase occurs for low temperature (T ~ 0.05) resulting in an arrested phase separation. Right 
panel: c = 4, I — 2 and s p = 2. The glass transition lines Td(p) and Tk(p) exhibit a reentrant shape (see inset), resulting in a 
glass-liquid-glass transition along isochores at densities around p ~ 0.62. 



The dynamical behavior is more involved [20(. The 
equilibration time diverges (as a power law) at the dy- 
namical transition pd — 0.5708 < px where the system 
is trapped in one of the exponentially many Gibbs states 
and falls out-of-equilibrium. When the pressure is fur- 
ther increased the system evolves within this state and 
the density can again be further increased [2l|, but does 
not follow the equilibrium curve. It is not yet known how 
to describe analytically this process on random graphs, 
and in FIG. Q] we thus display the results of Monte Carlo 
annealings at different rates. The slowest annealing we 
tried reached p f» 0.5731. We also computed the iso- 
complexity approximation curve for an infinitely slow an- 
nealing, where the dynamics is assumed to follow a curve 
defined by £(pDyn(/u)) = ^(Pd) [Hj]. This actually yields 
a lower bound on the infinite pressure density reached by 
this procedure. In our model, the density reached by an 
infinitely slow annealing, pj, is the analog of the Jam- 
ming point 23] and we thus obtain 0.5731 > pj > 0.5748. 
Note that, as discussed in 2l|, pj does not correspond to 
the densest possible packing /9rcp neither to the dynam- 
ical or Kauzmann density pd or px- In finite dimensional 
systems the dynamical transition is smeared out by acti- 
vation processes that restore the crgodicity and prevent 
a complete freezing of the system [14| . In that case, an 
ideal glass transition (if any) should thus arise only at 
the Kauzmann point pk ■ 

Phase diagrams - We now discuss the different 
phase diagrams obtained by varying the connectivity c, 
the maximum number of occupied neighbors i, and the 
entropic penalty s p . Three representative situations are 



depicted in FIG. O 

The left panel of FIG. [2] shows a typical phase diagram 
when the valency restriction I is small (but larger than 
2) and the dimension-like parameter c is large. In this 
situation a colloid-rich/colloid-poor phase coexistence is 
located in a narrow region of small temperatures T and 
densities p, in agreement with [9(]. The dynamical glass 
transition line Td then bends towards low value of p and 
continues in an almost horizontal direction. Under these 
conditions dynamical arrest phenomena and glassy be- 
havior should emerge at very small volume fractions. 
This state of matter is as close as possible to an ideal 
gel 0, which is in our model equivalent to a glass in 
the dynamical phase. In contrast, the Kauzmann glass 
transition line Tr- typically stays almost vertical in the 
T-p plane: a true equilibrium glass transition thus only 
takes place at high densities. In a finite dimensional sys- 
tem, where the dynamical transition is smeared out, the 
dynamics along isochores should follow approximatively 
an Arrhenius dependence, exactly as found in [1, Q. In- 
terestingly enough, the percolation line is located much 
above the dynamic arrest (gel) one. This gap is due to 
finite life time of bonds at finite temperature, this is com- 
mon for colloidal gelation. Nevertheless, the formation of 
a spanning cluster could affect the dynamical behavior of 
the system, as shown in [2J]. 

The central panel of FIG. [2] shows a typical phase dia- 
gram when the valency restriction I is large and c small. 
The coexistence region is much broader with respect to 
the previous case. And both the dynamical and the 
Kauzmann transition lines bend towards low densities. 
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An ideal glass state thus emerges for values of p much 
smaller than /Crop- The glass transition lines moreover 
intercept the spinodal so that a deep quench of the sys- 
tem at low enough temperatures (e.g., T 0.05) into 
the coexistence region (e.g., p ~ 0.2) results in a phase 
separation between a colloid-poor and a glassy colloid- 
rich region. Under these conditions the phase separation 
process is interrupted by the formation of an (attractive) 
glass in the dense phase. This arrested phase separation 
scenario is one of the possible routes to colloidal gelation 
and the formation of a low density gel state 

The right panel of FIG. [2] shows a case with the va- 
lency restricted to I = 2. In this case the occupied sites 
can create only chains and the phase separation region is 
absent, this was also observed in A more interesting 
feature shown in this diagram is the reentrant shape of 
the dynamical and the Kauzmann transition lines. The 
reentrant behavior can be observed when the short-range 
attractive potential is in competition with the entropic 
repulsion, s p , in agreement with At large temper- 
atures the interaction contributing to the formation of 
a glass is the entropic repulsion s p , at low temperature 
it is the attractive potential. Note that the reentrant 
behavior is usually observed together with the presence 
of the phase separation region, and we indeed observe 
the two effects in the left and central panel if we set the 
entropic repulsion s p to a positive value. Our example 
shows however that the reentrant behavior is not nec- 
essarily connected with the phase separation, nor with 
a presence of a glas s-glass transition, in agreement with 
earlier works [3. |25|. 

Conclusions - We have studied a lattice model of 
complex fluids, which can be solved exactly on random 
graphs using the cavity method [l5j |. Despite its intrin- 
sic simplicity, it captures many features that arise in real 
systems and it illustrates physical properties of arrested 
states observed in colloidal suspensions under different 
conditions. It reproduces the shape of a variety of phase 
diagrams, allowing for a qualitative comparison with ex- 
periments. In the future, it should be interesting to study 
the coarsening dynamics associated to the phase separa- 
tion between the colloid-poor and the colloid-rich phase 
and to simulate the dynamics in finite dimension, as in 
[lH l26l . [2^ 1 . It is also of interest to understand if a glass- 
glass transition can be found for this class of models. 

We thank A. Coniglio, J. Kurchan, O. Ri voire, T. 
Mora, F. Zamponi and E. Zaccarelli for discussions. 
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